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Abstract. - Solitons are studied in a model of a fiber Bragg grating (BG) whose local 
reflectivity is subjected to periodic modulation. The superlattice opens an infinite number of 
new bandgaps in the model's spectrum. Averaging and numerical continuation methods show 
that each gap gives rise to gap solitons (GSs), including asymmetric and double-humped ones, 
which are not present without the superlattice. Computation of stability eigenvalues and direct 
simulation reveal the existence of completely stable families of fundamental GSs filling the new 
gaps - also at negative frequencies, where the ordinary GSs are unstable. Moving stable GSs 
with positive and negative effective mass are found too. 



Bragg gratings (BGs) are distributed reflecting structures produced by periodic variation 
of the refractive index of an optical fiber or waveguide. Devices based on fiber gratings, such 
as dispersion compensators, sensors and filters, are widely used in optical systems [1]. Gap 
solitons (GSs) in fiber gratings arc supported (in the temporal domain) by the balance between 
BG-induced linear dispersion, which includes a bandgap in the system's spectrum, and the 
Kerr nonlincarity of the fiber material. Analytical solutions for BG solitons in the standard 
model of the fiber grating, based on coupled-mode equations for the right- and left-traveling 
waves, are well known [3]. Solitons with positive frequency w are stable, while those with 
lo < have an instability to perturbations that is characterized by a complex growth rate [4] . 

Solitons in fiber gratings have been created experimentally, with spatial and temporal 
widths of the order of a few mm and 50 ps, respectively [5]. Spatial GSs were observed 
in photorefractive media with an induced photonic lattice [6], and in waveguide arrays [7]. 
(Indeed, the new families of GSs reported in the present work may be realized in the spatial 
domain too, in addition to their straightforward implementation as temporal solitons in fiber 
gratings). Besides their relevance to optical media, GSs were also predicted [8] and created [9] 
in a Bose-Einstein condensate (BEC) trapped in a periodic potential. 
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An issue of particular technological importance is the development of methods for the 
control of BG solitons - in particular, using apodization [5, 10,11], i.e. gradual variation of 
the grating's reflectivity along the fiber. In an appropriately apodized BG, one can slow down 
solitons and, eventually, bring them to a halt [11]. Experimentally, it has been demonstrated 
that apodization helps to couple solitons into the grating [5]. Moreover, technologies are 
available that allow one to fabricate fiber gratings with periodic apodization, thus creating an 
effective superlattice built on top of the BG [12]. 

An asymptotic analysis of light propagation in such a superstructure grating was developed 
in Ref. [13]. It was shown that the superstructure gives rise to extra gaps in the system's 
spectrum ("Rowland ghost gaps"). Solitons in the gaps were sought by assuming that the 
soliton is a slowly varying envelope of the superlattice's Bloch function, which applies to the 
description of GSs near bandgap edges. Related problems were considered in Refs. [14], which 
treat BEC models with a doubly periodic optical lattice that opens up an additional narrow 
"mini-gap" , in which stable solitons may be found. 

The subject of this Letter is the investigation of GSs in harmonic (sinusoidal) supcrlattices 
created on top of the ordinary BG in fibers with Kerr nonlinearity. The analysis of the bandgap 
structure in this model constitutes, by itself, an interesting extension of the classical spectral 
theory for the Mathieu equation [15]. We will demonstrate that various families of GSs exist in 
the superlattice. Most importantly, in each newly opened gap we find a family of fundamental 
symmetric GSs, which fill the entire gap and are completely stable. Remarkably, they are 
stable not only at positive frequencies to, but also at 10 < 0, where the ordinary GSs are 
unstable [4]. 

The superlattice may be implemented via the creation of beatings in the optical interference 
pattern that burns the grating into the fiber's cladding. A model corresponding to this physical 
situation may be derived using standard coupled-mode equations for the amplitudes u and v 
of the right- and left-traveling electromagnetic waves. Inclusion of both Kerr nonlinearity and 
Bragg reflection terms leads to the following dimensionless equations: 

iu t + iu x + [1 — £ cos (kx)] v + [i cos (kx + 5 ) u + (\v\ 2 + \u\ 2 /2) u = 0, , , 

ivt — iv x + [1 — e cos (kx)] u + fx cos (kx + S) v + (\u\ 2 + \v\ 2 /2) v = 0. 

Here e is the amplitude of the periodic modulation of the BG strength, while the nonlinearity 
coefficient and average reflectivity are normalized, as usual, to be 1. Accordingly, k measures 
the ratio of the modulation period L to the BG reflection length I (usually, L < 1 cm [12, 13] 
and I ~ 1 mm, while the underlying BG period is ~ 1 fira, whereas the total length of the 
grating may be up to 1 m). We have included additional terms in Eqs. Q with amplitude /i 
and phase shift S to describe another possible control mechanism for optical pulses in BGs [1,2], 
namely a periodic chirp of the BG, i.e. a local variation of the grating's period. 

In this work we focus on the reflectivity superlattice, setting ji = 0; the general situation 
including the periodic chirp will be considered elsewhere. Thus e and k are the parameters of 
the model, which we assume to be positive without loss of generality. According to the above 
physical considerations, typically e is small, but k is not. Note that we should only expect a 
quiescent GS to be stable in this kind of nonuniform BG if the position of the soliton's center 
is located at a local minimum of the reflectivity [11], i.e. at x = 2wn/k for some integer n. 
Hence we shall restrict our numerical search for GSs to just such a case. Lastly, note that 
Eqs. JjQ| conserve two dynamical invariants, the Hamiltonian and norm (usually called energy 
in fiber optics), E = f_™ (\u(x)\ 2 + \v(x)\ 2 ) dx. 

We shall seek stationary GS in the form {u(x, t), v(x, t)} = cxp(— iu>t) {U(x),V(x)}, sub- 
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Fig. 1 - The spectrum of Eq. Q with k = 1 in the (w,e) plane. Shaded and white areas are, 
respectively, bands and gaps. Digits label the gaps mentioned in the text. Markers correspond to the 
solutions plotted in Fig. [5] 



stitution of which into Eqs. 0) yields 

uU + iU' + [1 - e cos(kx))V + (\V\ 2 + \U\ 2 / '2) U = 0, (] 
ujV -iV' + [l-e cos(kx)} U + (\U\ 2 + \V\ 2 /2) V = 0, ( ' 

where prime stands for d/dx. First, it is necessary to find bandgaps in the system's spectrum. 
In the unperturbed problem (e = 0), the dispersion relation for linear waves with a propagation 
constant g is w 2 = q 2 + 1, thus producing the well-known bandgap, — 1 < uj < 1. The 
spectrum for e > should be found from the linearization of Eqs. J5J). As in the Mathieu 
equation, gaps in these equations emerge due to the parametric resonance caused by the 
cosine modulation. Straightforward considerations show that new bandgaps open up at points 
lo = uj m = sgn(m)i/l + (mfc) 2 /4, with m = ±1, ±2, . . . (the extra modulation terms ~ ji in 
Eqs. Q] ) open gaps at the same positions). We will designate these gaps as m.^. Using 
perturbation theory, one can find their boundaries for small e and demonstrate that their 
widths scale as e' m l. In particular, gaps 1^ are found in the region |x| < 1, x = 2 (w — uj±i) /e, 
to leading-order in e. 

Equations @ allow the invariant reduction V = —U*, which reduces them to a single 
equation, 

loU + iU' - [1 - e cos(fcr)] U* + (3/2)\U\ 2 U = (3) 

(the reduction with V = U* leads to an equation with the opposite sign in front of U* , which 
can be cast back in the form of Eq. (|3J by the substitution U = iU; neither reduction is valid 
for moving solitons, see below). The linearization of Eq. is sufficient for the analysis of the 
spectrum. Bandgap regions in the (u>, e)-parameter plane were computed using the software 
package AUTO with driver HomMap [16]. We detected points at which the linearization of 
the Poincare map [17] around the origin (U = 0) has eigenvalues ±1, and then continued 
such points to identify the bandgap boundaries. The results are displayed in Fig. ^ which 
shows the first five gaps. In this and subsequent examples, we set k = 1, as this value 
adequately represents the generic situation and is physically meaningful for the application to 
fiber gratings. 

In the central gap, which we designate 0, standard calculations using Mclnikov's method 
[17] reveal that the GS solutions, which are known in an exact form for e = [3], extend to 
e > 0. They were continued numerically up to e = 1, and were found to completely fill the 
central gap. That is, a GS exists for each {w,£}-value belonging to the gap. We have also 
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computed eigenvalues that determine stability of these solutions against small perturbations, 
and found that only GSs with ui > arc stable. That is, the border between stable and 
unstable GSs in the central gap for e > is found to remain, up to numerical accuracy, at 
u) = as it is for e — [4] . 

Inside gaps 1 ± , one can use the method of averaging [17] to demonstrate the existence of 
GSs for small e. To this end, we represent a solution to Eq. J3J, U(x) = a(x) + ib(x), as 



= V^exp 



-(w±i + l) 

LU ±1 - 1 



Elk 

e/(2(w±i + i)) 



(4) 



With constant amplitudes S and O, we recover a solution to the linearized equation with 
e = and u) = uj±. The averaging method supposes that E and O are functions of a slow 
coordinate, z = xj [2k(uj±i + 1)], which leads to equations 



±( e\_( «(i - X )Q - p (e 2 + e 2 ) e 

dz 1 6 J ~ I a(l + X )E + (3 ~ 2 + e 2 5 



(5) 
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where a = e [k 2 + 4(w ±1 + l) 2 ], and/3 = 3e[3k 4 + 8(oj ±1 + l) 2 k 2 + 48(w±i + l) 4 ]/[8fc(w±i + l)] 

These equations conserve their Hamiltonian, H = a [(1 + x)S 2 — (1 — x)© 2 ] +(A/2) (^ 2 + 6 2 
As gaps 1 ± exist for |%| < 1, the coefficients 1=FX m Eqs. (JSJ and in 7f are positive. Therefore, 
the origin (S, 0) = (0, 0) is a saddle in Eqs. 10), and a pair of homoclinic orbits to this saddle 
can be found in the exact form 

(E±(t),G±(t)) = ±2 J— sin6»sin- — - \ ns ' -, 6 

v w w; V V 3/3 2 cosh(2atsin6») +cos6> w 

! cosh(atsin^) \ 

C ° b 2 cosh(2at sin 6) + cos 6 ) ' 

with x = cos0. On application of the transformation ijlf. these orbits correspond to solitons 
in the full system In particular, there are solutions with a(x) = KeU(x) odd and b(x) = 
ImU(x) even. Thus, GSs exist in the entire gaps ^ for sufficiently small e. 

A similar analysis can be performed for higher bandgaps, 2 ± , 3 ± , etc., using a higher- 
order averaging method [18]. The respective analytical computations show that each averaged 
system again generates solitons. They have either even real and odd imaginary parts or vice 
versa, depending on whether the gap's number m is odd or even. 

We employed the AUTO driver HomMap [16] to continue numerical soliton solutions of 
the full system @, varying lu and s . In so doing, a multitude of symmetric and asymmetric 
families of GS solutions, including higher-order ones (bound states) were found in each new 
gap. Here, we display numerical results for fundamental solitons, as their bound states are 
likely to be unstable. It was found that, as predicted above, each gap 1 ± , 2^, etc. is completely 
filled with a single family of symmetric GSs. The families in gaps 1 ± are represented by 
solitons displayed in Figs. Gl^) and (b ± ), for small and larger e, respectively, while Figs. 
I2c ± ) displays the fundamental GSs in gaps 2 ± . Note that stable fundamental solitons may 
feature a double-humped structure in terms of |E/(x)| at relatively large e, which is the case for 
the soliton in Fig.|2Ib + ), sec Fig. Efb) below (however, the |t/(x)| profile of the GS displayed, 
for the same e but opposite u>, in panel I3b _ ), remains single-humped, see Fig. Ota)). The 
homoclinic orbits of averaged system (J5J, also shown in a) and c), provide a good match to 
the envelopes of the numerical solutions. 

Recall that the ordinary GSs in gap with e = cannot be asymmetric, and they do 
not form bound states either [2] . To explain where the asymmetric and higher-order GSs for 
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Fig. 2 - Stable fundamental soliton solutions of Eq. © in gaps 1* and 2 . Gray lines in rows a) and 
c) show solutions of averaged equations @, and their counterparts for gaps 2 . 



£ > come from, we note that the new symmetric fundamental GSs, found above in gaps 
1^, correspond to transversal intersections of stable and unstable manifolds in the Poincare 
map associated with Eq. (J3 ) . Such intersections naturally occur in pairs, one representing a 
symmetric soliton and the other one its asymmetric counterpart. Moreover, the transversality 
of the intersection implies the presence of a Smale horseshoe [17], which, in turn, ensures 
the existence of infinitely many higher-order GSs, that may be realized as bound states of 
symmetric or asymmetric fundamental solitons. 

Stability of the GSs was tested by direct simulations of Eqs. JU, and verified through 
computation of the eigenvalues of the equations linearized around a GS, which govern the 
growth rate of small perturbations. Using this approach, it was found that the entire families 
of the symmetric fundamental solitons that fill the gaps 1 + and l - are stable. The significance 
of this result is that we can have stable GSs with to < 0, something, that is not possible in 
the ordinary BG model without the periodic modulation terms. This finding is illustrated in 
Fig-Efa), taking as an example the soliton with u = —1.118 depicted in Fig.[2Ib~). 

Regarding asymmetric solitons in gaps with m > 1, it was found that some of them 
are stable and some are not, the unstable ones being completely destroyed by growing pertur- 
bations (not shown here). 
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Fig. 3 - Soliton dynamics in gaps 1 T . (a) Evolution of a perturbed symmetric soliton shown in Fig. 
Hfb - ) (e = 0.4, uj = —1.118; for these e and uj, an asymmetric soliton exists too, but it is unstable), 
(b) Splitting of the stable double-humped soliton shown in Fig.|5Jb + ) into moving and quiescent ones, 
after it was suddenly multiplied by exp (ipx), with p = 1. 

An important issue for practical application is the possibility of finding moving GSs. (In 
experiments on ordinary fiber gratings, only moving solitons have so far been observed [5]). 
In the presence of the superlattice, soliton mobility is a nontrivial problem, because the GS 
has to move in a periodically nonuniform medium. (Recall, though, that solitons belonging 
to gaps m ± with m > 1 do not exist at all without the superlattice). We have used numerical 
simulation to test for mobility of the stable GSs we have so-far found. Specifically, a quiescent 
stable soliton was multiplied by exp(ipx), with p > 0, which implies a sudden application of 
a "shove" to the soliton, giving it momentum P = i (uu* + vv*) dx + c.c. = 2pE, where 
E is the soliton's energy defined above. It was observed that GSs belonging to central gap 
are readily set in stable motion by the shove. This accords with the situation for e = 0, 
where exact solutions for moving GSs exist with any velocity c from — 1 < c < +1 [3]. In 
our simulations, we found that no soliton in any gap could be made to move with velocity 
exceeding 1 (a very strong shove does not make the GS "superluminal" , but simply destroys 



For the stable GS belonging to gap 1 + , it was found that the application of the shove 
caused the soliton to split into a quiescent part and a moving part; see Fig.^b). A remarkable 
observation from these results is that the moving soliton (which retains ~ 2/3 of the initial 
energy) has a negative average velocity, c « —0.9, hence its effective mass M ~ P/c is negative 
too. (Simulations of the moving GSs belonging to gap have positive mass, as is the case for 
e = [11]). We remark that the mass of a GS may be negative too in one- and two-dimensional 
models of BEC in an ordinary optical lattice, and that this gives rise to nontrivial effects such 
as stable confinement of solitons in an anti-trapping external potential [19]. Finally, the 
shove applied to a stable GS belonging to gap 1 was found to split into three solitons, one 
remaining quiescent while the other two move off in opposite directions (not shown here). 
Detailed results for moving solitons will be reported elsewhere. 

In conclusion, we have investigated a model for a Bragg grating optical fibre with a super- 
lattice written on top. First, we identified a system of new bandgaps in the fiber's spectrum, 
using an extension of the bandgap theory for the Mathieu equation. Combining averaging 
methods and numerical continuation, we have found that each new gap is completely filled 
with fundamental symmetric solitons, which at some parameter values may have a double- 
humped shape. In addition, new types of gap solitons (GSs) were found that do not exist 
without the superlattice, such as asymmetric GSs and bound states of fundamental GSs. An 
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important finding is that the entire families of fundamental GSs in the new bandgaps are sta- 
ble, including negative frequencies, where ordinary GSs are unstable. Stable moving solitons 
were found too, including ones with a negative mass. Finally, it is pertinent to point out that 
creation of the newly predicted GSs ought to be perfectly feasible using presently available 
experimental techniques. In particular, in a weak superlattice, with say, e ~ 0.1, an estimate 
shows that the solitons in the new bandgaps l^ 1 can be observed if the fiber grating, on top 
of which the superlattice is to be imposed, is of length ~ 10 cm. 
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